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ABSTRACT 

The average flow model of Patir and Cheng [1,2] for obtaining an aver- 
age Reynolds equation in the presence of two-dimensional surface roughness 
is extended and generalized. Expectation values of the flow factors appear- 
ing in the formalism are calculated by means of a perturbation expansion of 
the pressure in a nominal parallel film. Terms in the series are evaluated 
using the unperturbed Green function which permits ensemble averaging to be 
performed directly on the solution. Calculations are carried to second 
order which involves only two-point correlation functions of the two rough 
surfaces. Perturbation results agree well with results of the earlier nu- 
merical simulation until surface contact becomes important when both 
approaches are inadequate. The theory displays the dependence of the flow 
factors on the roughness parameters in simple closed form, leading to 
improved understanding of the average flow method. 

INTRODUCTION 

In the theoretical study of the effects of surface roughness on the 
performance of lubricant films, the objective is to develop a form of the 
Reynolds equation from which the expected pressure at each point in the 
fluid may be determined. While limitations on the ability of the Reynolds 
equation to describe flow between rough boundaries continues to be a topic 
for much discussion [3], the wide variety of practical surface textures 
falling within the generally accepted domain of validity nonetheless justi- 
fies this endeavor. Among the basic equations of mathematical physics, the 
Reynolds equation is unusual in that boundary conditions on the flow deter- 
mine the form of the equation itself, rather than serving only to define a 
particular solution of a single equation. These boundary conditions, giving 
the shape and velocities of the surfaces, enter as a result of the partial 
integration of the Navier-Stokes equation across the thin dimension of the 
film. In the roughness problem, where the shape of the bounding surfaces is 
a stochastic process, it is this feature which offers the possibility of 
constructing the required ensemble-averaged form. It is not the purpose of 
this paper to discuss the many previous attempts to develop such equations 
of which several critical reviews already exist [4]. However, as noted by 
Elrod [3], most treatments of the Reynolds equation have not encompassed the 
general situation of two-dimensional roughness, with consequent limitations 
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on their applicability. Prior to Elrod's own work, one of the first attacks 
was made by Patir and Cheng [1,2]^ who proposed an ensemble-averaged 
Reynolds equation in which roughness effects were built into a number of 
flow factors. Since their pioneering study, many authors [5-7] have made 
use of the results either directly or for comparative purposes and it is 
timely to establish the PC approach on a more analytic footing. The present 
work constitutes both a critical analysis and an extension of their average 
flow method, concluding with a discussion of some of the conceptual problems 
encountered in the partial lubrication regime and an indication of how such 
problems are being tackled. 

In the PC formulation the actual flow between rough surfaces is equated 
to an averaged flow between nominally smooth surfaces, while parameters de- 
scribing the roughness are included in the Reynolds equation through the 
flow factors. The problem is thus reduced to an evaluation of these fac- 
tors, accomplished by considering flow between nominal plane parallel bound- 
aries. At this stage, PC adopt the direct approach of numerical simulation 
of the rough boundaries together with numerical solution of this model flow 
problem. Apart from the inherent difficulty associated with extracting 
reliable expectation values (ensemble averages) of the flow factors found 
for each sample roughness, this approach is also unsatisfactory in that the 
functional dependence of the factors on the roughness parameters can only be 
inferred from numerical experiments aid expressed in fitted form. 

We overcome these difficulties here by developing formal solutions of 
the model Reynolds equation which may be explicitly ensemble-averaged. This 
automatically introduces relevant statistical properties of the film bound- 
aries. The analysis not only agrees quantitatively with the simulations of 
PC but also places on a firm theoretical footing several conjectures made in 
that work on the basis of the numerical experimentation. 

The formal solution is written in terms of the Green function for ideal 
smooth flow, enabling the pressure to be expressed in a form containing the 
random height of each surface. The flow factors themselves, involving pro- 
ducts of film thickness and pressure derivatives, may thus be directly aver- 
aged with no appeal to an ergodic assumptionc A perturbation series is 
developed for the pressure in which the expansion parameter effectively 
'turns on' the roughness, with zeroth order corresponding to smooth flow. 

One important variable, H, governing the flow is the ratio of the nominal 
film thickness to the standard deviation of the combined roughness. For 
H ^ 3 areas of contact near surface summits become progressively more sig- 
nificant both in blocking flow and in altering the statistical properties of 
the surfaces. In this partial lubrication regime, our results begin to 
deviate from those of the PC model, which in contrast to the present work 
makes some allowance for contact. 

NOMENCLATURE 


g-j smooth pressure gradient 

G Green function, solution of Eq. (10) 

h film thickness = Z 2 - zj 

H non-dimensional nominal film thickness = h* / a 


1. These two papers will be referred to together as "PC". 
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Cartesian tensor indices 
film pressure 

tensor determining pressure flow, Eq. (16) 
reduced film pressure = (1 - e““P)/a 
volume flow 

right-hand side of Eq. (Al) for q^ 
tensor determining shear flow, Eq. (16) 
time coordinate 

surface translational velocities 
mean surface velocity = (ug + 
slip velocity = U 2 - u^ 

Cartesian space coordinates 
heights of bounding surfaces 1 and 2 
mean surface height = {z\ + Z 2)/2 

viscosity pressure index 

correlation function parameters, Eq. (19) 

roughness anisotropy index (Peklenik number) = (s/a)^/2 

Dirac delta function 

Kronecker delta = unit matrix 

roughness heights of bounding surfaces 

perturbation expansion parameter 

non-dimensional sum and difference roughness heights = 
(«2 * fil)/h* 

coefficient of viscosity 

generalized two-point correlation function 

two-point correlation function of (62 - s^) 



a 

rms value of (62 - «i) 

ol,2 

rms values of «x,2 

^x»i^y»^s 

x,y pressure and shear flow factors given by PC 

*^ij» ^ij 

pressure and shear flow tensors 

♦s 

single surface shear flow factor 


Special Symbols 
★ 

expectation value operator for stochastic quantity 


FORMALISM 


The class of flow problems considered by PC is illustrated in Fig. 1 
which shows a fluid film between two rough boundaries, 21 and zg, 

moving in the (x,y) plane with arbitrary translational velocity, ui and 
U2 measured relative to a point in the fluid remote from the film region. 
This specification of surface velocity eliminates rigid rotation, although 
in cases where the nominal surfaces^ zx* and Z2* possess rotation axes fixed 
relative to one another, a straightforward extension of the formalism 
applies. The actual surface height zi 2 is considered to be simply the 
sum of the specified nominal height and’a random roughness height 6x 2* 
so that difficulties in choosing a wavelength cutoff for roughness arS 
avoided. The random variable 6^ 2 a stationary stochastic process 
distributed with standard deviatibn ai 2 about zero mean. We need also 
the standard deviation o of the rougnftess combination (62 - fix)» 
which in cases where cross-correlation is absent will be just 
2 2 1/2 

(oj + 02) . For convenience we replace the boundary variables with the 

following combinations: 


film thickness h = Z2 - zx; mean surface height Z = (Z2 zx)/ 2 ; 

slip velocity 'y = U2 - ux; mean surface velocity 0 = (u2 ui)/2. 

Making the usual thin film assumptions, the volume flow per unit width through 
the film is 




a.p + h U. 


( 1 ) 


2 . We use the asterisk to denote nominal or expectation values (ensemble 
averages) except where several factors are included when the angle bracket 
<> notation is clearer. 
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where p is the film pressure and y the shear viscosity coefficient of 
the lubricant.^ For incompressible flow the continuity equation 
3iQi = -8h/3t is applied in the usual way to derive the Reynolds equa- 
tion for this problem 
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-Vi3iZ 


( 2 ) 


where the Euler transform has been used to interchange time and space deri- 
vatives of h. 

From Eq. (2) it is clear that p is a stochastic variable correlated 
with h so that care is required in calculating average flow from Eq. (1). 
The goal Is to write in terms of h* and 3^p* so that applying 
continuity will yield a^Reynolds equation in which only averaged quantities 
appear. There is some latitude in the way this step is taken and we follow 
FC as closely as possible. The Couette term presents no problem and is 
handled by replacing h with h* so that changes in the form of 
Q* arise only from the pressure term. We note from Eq. (2) that even in a 
nominal parallel film there will be a term in in the pressure gradi- 
ent, while likewise, even if V^- vanishes the cofactor of the pressure 
gradient is modified. We separate these two effects on the pressure and 
write the flow as 


Q 


i 


*3 

-h 


.3 .p* + hu. - 


7 .V . 

2 '"ij j 


(3) 


This equation serves to define the pressure and shear flow factor 

tensors and respectively. In general these tensors are symmetrical 

but not diagonal so that, for example, an x pressure gradient can produce 
a y flow under certain conditions of roughness anisotropy. While Eq. (3) 
is exact, its utility rests on the hypothesis that the relative variances of 
the flow factors are small compared to those of h and p, in which case 

the Reynolds equation derived from Qt instead of Q- will represent roughness 
effects satisfactorily. The basis for this hypothesis lies in Eq. (3) which 
shows that Q-j and the ^ factors have the same relative fluctuations. 

Since continuity determines flow globally rather than locally, its fluctua- 
tions are indeed smaller than those of h or p. It i| incorrect, however, 
to assume flow fluctuation is exactly ze^o when only h is specified and 
a distinction remains between i and 0 . The PC factors correspond to 
the latter which by comparing ensemble-averaged Eqs. (1) and (3) must satisfy 





3 .p + . V . 


(4) 


3. Henceforth subscripts i and j refer exclusively to Cartesian compo- 
nents with (xi,X 2 ) denoting (x,y). Summation on repeated (i,j) sub- 
scripts is implied. 


DRfGINAl. IS 

OF POOR QUALITY 
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ORlGiNAi PAQjff fg 

Finally, using Eq. (4) in ensemble-averaged Eq. (2) yields POOR QUALITY 


^\^ij W “ ■* ''j®i^^*'^ij ■*■ 7 

where is the unit matrix (Kronecker delta). This average Reynolds 

equation is essentially that proposed by PC apart from the appearance of 
the off-diagonal flow factor components, allowing Eq. (5) somewhat greater 
generality. In the following section where we calculate the complete flow 
factor tensors, we compare components with those of PC according to the 


(5) 




s* 


the correspondence = ^x’^22 ^ ^11 ^ ^s * components 

do not appear in PC. From Eq. (4) we see that + 6 . . and 0^* + 0 in 

* w V ♦.* 

the smooth limit so that 0;^ y represents fractional increase in (x,y) 
flow in absence of slip, while the other c'.nponents represent various modes 
of entrained flow not present in the smooth case. 


CALCULATION OF FLOW FACTORS 


The flow factors are calculated from Eq. (4) by first solving the exact 
Reynolds equation (2) for p and then taking expectation values. We apply 
Eq. (2)^to the flow between nominally plane parallel surfaces with separa- 
tion h and develop the solutjon in perturbation series to obtain the 
flow factors as functions of h" and certain roughness parameters. In 
contrast, PC solve Eq. (2) numerically using a populatioi. of carefully 
generated rough surfaces. Their ensemble-averaging is then executed in two 
stages, first invoking the ergodic hypothesis to take a spatial average of 
the flow across the model bearing area for each sample, followed by a popu- 
lation average. Elrod [3] also requires the ergodic hypothesis in his per- 
turbation treatment. 

Before proceeding, one further simplification is possible if the 
pressure dependence of the viscosity is restricted to the exponential form 

u = reduced pressure q = {1 - e”“^)/a satisfies 

-1 -1 

a.q = p s^p so that p and p may be replaced by p^ and q. 

If the surfaces zi 2 represented by = - h*/2 + fii 
and Z 2 = h*/2 + 62 , then the variables h and Z become: 

h = h* + 62 " = h*(l+x“) and Z = (62 <si)/2 = fi*x'*‘/2 

which defines the non-dimensional sum and difference roughness coordinates 
= ('52'*' iSi)/h*. In the spirit of perturbation theory, we introduce 
an attenuation factor e, 0 < e _< 1, such that roughness effects are intro- 
duced by letting e >1. Using e as the formal expansion parameter we 

CD 

write q = 1C which, upon extracting the coefficient of converts 

n=0 " 

Eq. (2) into a set of equations for qp. Because of the factor h^, the 
general form occurs for n ^ 3 but here we are concerned mainly with special 
forms for n = 0 and 1. These are 
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( 6 ) 


n = 0 , 3^3^qQ = 0 

n = 1 . 8,3,qi = -3a,{x-s,q ) - (7) 

n* 

Up to terms in the expression required in Eq. (4) becomes 
l,3 *3 ( 

+ e^[3(\")^ 3jq^ + 3x"Sjqj + s.q^] + 0 (e^)| (8) 

Comparison of the expectation value of Eq. (8) with Eq. (4) will then yield 
the t tensors correct to second order in roughness. In the Appendix it is 
shown that <3^q^> = 3^q* = 0 for n > 0 so that roughness leaves the 

pressure gradient unchanged from its smooth value. Moreover from Eq. (6) 
this may be assigned the constant value, g-j, equal to the mean applied 
gradient. With this simplification, Eq. (8) may be averaged to give 

3^-^= T^jgi 3e^[H"^g. + <X~3.q^>] + 0(e^)| (9) 



Since q 2 no longer appears, to this 
tion is <x”3-jqi> which contains all 
roughness. 

The solution for qj is written 
of the smooth problem, which from Eq. 

»,3^G(Xj;xj) = 6(Xj-x]) 


order the only term requiring atten- 
the interesting effects of 

in terms of the Green function G 
(6) satisfies 

( 10 ) 


where 6 is the two-dimensional Dirac delta function. Ignoring the comple- 
mentary function, the particular integral of Eq. (10) is well known to be 

G(Xi;x|) s.n [(xpxf)(x. -x|)] (11) 

which has the important property of depending only on the difference of its 

arguments. Using the notation si- to denote 3/3x1-, the solution of Eq. (7) 
for q^ becomes ^ ^ 

q.(x) = -3g^jyG(x-x^)a|x~(x^)dx^ - ^^o \/.JJ G(x-x^)3|x^(x^)dx^ (12) 

where subscripts have been omitted from spatial variables. From this the 
expectation value appearing on the right of Eq. (9) is found to be 
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<x = 


3gjf|3^G(x-x^)3^x”(x)x“(xV(^x^ OR/ginaL PAGE IS 

Op POOR QUALiry 


- 


vjj3^G(x-x^)8j<x”’(x)x^(x^)> dx^ 


(13) 


This shows that all the required statistical information lies in the two 
two-point correlation functions (CF) defined here. The first CF, p^, 
is that for the combined roughness variable (62 - which in conven- 
tional form becomes 

<x"{x)x"(x^)> = H"^p^(x-x^) (14) 


using the definition of = h*3<x-(x)x"(x)>. It; many cases of 
interest, the cross-correlation between the two surfaces < 6 i(x) 62 (x 4 > 

2 2 2 

will be zero so that a = 

2 2 

°2 ^2 2 autocorrelation functions of the 

T" ~T- 

o a 

two surfaces. The second CF, pq , is related to a difference of the 
ACF's of the two surfaces and may be written 

<x~'(x)x^(x^)> = Pp(x-x^) . (15) 


By symmetry any cross-correlation vanishes from pn which thus becomes 
2 2 

^2 p« - ^1 f-., . In showing the dependence of these CF's on the argument 
~T~ ^ ^ 

a a 

difference, the stationarity of the roughness process has been invoked. 

It is useful to assign names Pt- and S-jj to the tensors defined 
by the integrals in Eq. (13) according to their association with pressure 
and shear flow, thus 



in which the difference argument no longer need appear. 
(13) and (16) and comparing with Eq. (4) then leads to 


(16) 

Combining Eqs. (9), 



ij 






(17) 




(18) 


This completes the formal part of the calculation. |he <f) tensors in 
Eqs. (17) and (18) lead via Eq. (3) to an averaged flow Qi structurally 
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similar to that derived by Elrod. The essential difference lies in the 
direct appearance of the CF's in Eq. (16) by contrast to their Fourier 
transforms (power spectral densities) in that earlier study. In other 
respects the Green function technique complements the Fourier transform 
method as a means of solving the Reynolds e'ljuation. 

RESULTS 

For numerical comparisons it is now necessary to construct some model 
CF's to insert in Eq. (16). As we have seen in connection with Eqs. (14) 
and (15)» pc and pp in the absence of cross-correlation are just a 
weighted sum and difference respectively of the individual ACF's of each 

surface and the factors are correspondingly expressed in sum or differ- 

ence form. Our discussion of the general form for the CF applies to any one 
of these four functions. Even where cross-correlation is involved we assume 
the form of the CF still falls within the general description. Hence, appro- 
priate combination of results for a single CF yields results for any given 

pair of surfaces. The flow factors in any case are not expected to be sen- 

sitive to the precise CF form provided they refer to physically similar sur- 
face textures. 

Our calculations correspond to textures having no long-range order, for 
which the CF necessarily decays to zero as the correlation distance 
increases in any direction from the central summit, p = 1, at the origin. 

In general, the contours of p possess only the symmetry of inversion 
through the origin but if any preferred direction exists as a result of the 
finishing process, the existence of two perpendicular planes of reflection 
symmetry for the CF is assured. In such a situation which includes isotropy 
as the special case where all planes are reflection planes, it is possible 
to define Cartesian 'roughness axes' and since no restriction has been 
placed here on the boundary velocities or pressure gradients, no generality 
is sacrificed by choosing (x,y) to coincide with these axes. This corre- 
sponds formally to a coordinate rotation which diagonalizes the flow factor 
tensors, the simplest representation from a computational viewpoint. When 
roughness axes exist, the Peklenik y factor [8] is well-defined, being the 
ratio of the x and y intercepts of the p = 1/2 contour. While any y 
value still embraces an infinite variety of surface textures, it is a useful 
characterization adopted by a number jf authors, including PC, to describe a 
surface whose asperity contours typically have an x dimension y times 
their y dimension. In many cases these will be the maximum and minimum 
dimensions. The special values 1, 0 and «> correspond respectively to 
isotropic roughness and one-dimensional ridges parallel to y or x. 

Conforming to these most general specification, PC chose a CF decaying 
linearly with distance and having rectangular contours. This form is con- 
venient fo, their numerical simulation. We have performed the integration 
of Eq. (16) for this choice and find that the results differ unimportantly 
from a Gaussian form for p which has the advantage of a remarkably simple 
final expression for the flow factor components. Thus, for purposes of 
illustration we adopt 


p = exp - [a(x-x^) + B(y-y^) ] (19) 

for which the contours are elliptical with y = V s/a . With this form for 
p, the range of integration in Eq. (16) may be extended to infinity with in- 
significant loss of accuracy, as discussed for example by Bush and Gibson [7] 
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in a Green function analysis of the inclhied slider bearing. Substituting 
Eq. (19) into Eq. (16) for pr then leads to the result Pn (v) « -(y+I)”^ » 
while P 22 (y) ° Pn (1 /y)* Thus, in PC notation 

<^j^(h,y) ° 1 + 3^ original page is 

OF POOR QUALITY 

^y(H,Y) = fij^(H,l/Y) . (21) 

Similarly if p denotes Pj^ or pg , we obtain 

0 ^ ^2 

<^s(H,Yi,Y2) “ (*^»Y^) 2 — 

a o' 

where 

»s(H,-r)=|^ (23) 

is the shear flow factor for a single rough surface. These closed expres- 
sions may now be compared directly with the PC numerical results. 

DISCUSSION 

In Fig. 2 are shown some plots of Eqs. (20) and (23) together with the 
numerical results of PC. As we should expect, there is insignificant dis- 
agreement between the two as H grows large and effects of roughness on 
flow diminish. Differences being to appear for H<3 in the partial lubri- 
cation regime. Since contact effects are not included in the present work, 
results are plotted only for H>2. 

This calculation has several interesting features. The asymptotic 
agreement of numerical simulation with the perturbation approach to solution 
of the Reynolds equation indicates that both techniques are accurate in the 
absence at least of significant degree of contact. In the two special one- 
dimensional cases, for which independent estimates of the flow factors are 
available [9], the perturbation method again agrees asymptotically. It thus 
supplies a useful analytic expression connecting these two extremes. The 
most striking feature of the pressure flow factors exhibited by PC is the 
crossover from monotonic increasing to decreasing behaviour occurring for a 
Y value between 1 and 3. Eq. (20) indeed suggests the transition occurs 
exactly at y 2 for which '’s constant. Although this exact con- 
stancy depends on the choice of model CF (with the linear CF used by PC, 
dx is found to increase very slightly for y = 2), this qualitative 
agreement is nonetheless most encouraging. It sheds some light on the 
physical basis of this curious transition from reduced to enhanced flow as 
the longitudinal contour dimension parallel to the pressure gradient) 

exceeds twice the transverse. Eq. (20) shows that transverse roughness is 
twice as effective in reducing flow as longitudinal roughness is in 
enhancing it. Thus if cross-correlation is zero, a combination of trans- 
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verse furrows on one surface with longitudinal furrows on the other has no 
resultant i?ffect on flow = 1) » provided the variance of the longitu- 
dinal surface is just twice the transverse variance. 

In general, the distinction between our flow factors and those of PC 
is that in each case ours predict a relatively larger roughness effect. 

Thus, our 4's factor inevitably must be larger at small H since contact 
does not cut off the flow, and likewise the enhancement of pressure flow for 
longitudinal types of roughness will be overestimated. It is not so easy to 
interpret the present overestimation, at least relative to PC , of the 
reduction in flow produced by transverse texture. It is, however, clear 
that as Y decreases, fluctuations in calculated by simulation 
become more important and may make such techniques unreliable. 

The analysis presented here also shows clearly how the flow factors are 
determined by the combined stochastic properties of the two surfaces. The 
appearance of two distinct CF's, pc and pg in pressure and shear 
flow respectively, is an important feature not discussed by PC. Whereas 
Eq. (22) here follows directly from Eq. (16), PC are only able to infer it 
from numerical experimentation. Similarly we can demonstrate another con- 
jecture that so long as the y factor is uniquely definable, the precise 
form of the CF does not strongly influence the numerical results. As 
already indicated, integration of Eq. (16) for the linear ACF used by PC 
alters our results by much less than the differences apparent in Fig. 2, 
while using yet another CF , Elrod's results agree closely with ours [3]. 
The idea receives further support in the results of Bush and Gibson [7J, 
while the power of the form of the CF to predict tribological parameters 
in general has been emphasized in the work of Whitehouse and Phillips [10]. 

OUTLOOK 

At this point the chief advantage of the Green function perturbation 
approach to the flow problem in general becomes apparent, namely the 
explicit way the statistical properties of the surfaces appear. Continua- 
tion of the perturbation series for pressure to arbitrary order n leads, 
as shown in the Appendix, to the involvement of all t-point CF's up to 
= n. Insofar as such CF's are experimentally measurable, their informa- 
tion could be included directly in the calculation of flow factors which 
goes some way to answering T(&nder's objection [2,5] that the parameters a 
and Y insufficiently characterize roughness. Even to second order there 
are no restrictions beyond the general ones discussed by the form given to 
the CF's. The formalism, however, is more difficult to apply to surfaces 
with long-range order of the type constructed by Tjinder which have no 
unique y value. 

The chief cause of breakdown of the present approach, however, is not 
truncation of the perturbation series but the changes in statistical proper- 
ties resulting from contact growth at small H. Contact not only alters the 
surface ACF's but also introduces cross-correlation. At higher degrees of 
contact the most important effect for flow is the shape and connectivity of 
the contact spots. Meeting the required condition of zero flow normal to 
any contact boundary, together with zero flow within pools of fluid trapped 
inside a closed boundary, involves the topology of the entire interface. 

The PC method which consists of setting local flow zero whenever the sim- 
ulated heights indicate surface overlap makes some allowance for contact but 
grows progressively more inadequate as H decreases. It is likely that 
for H 2 a percolation description of flow is more appropriate. Percola- 
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tion theory, for example, demonstrates that when the fraction of contact 
area In a random isotropic (y = 1) case reaches 44 percent, there no longer 
exist any open paths for flow [11]. On a Gaussian surface this occurs for 
H ~ 0.15. Rather surprisingly then the flow factors should vanish altogether 
at this point, even though the expectation of the film thicknes*’ (ji;ean gap) 
is nonzero. 

Such an approach to the partial lubrication regime emphasizes the ran- 
dom nature of the contact geometry. An approach described by Shukla [12] 
with somewhat the same viewpoint treats the interface as a random network, 
presenting an effective permeability to the lubricant, the apparent viscos- 
ity of which depends on the degree of contact. Such methods appear well 
suited to extending our understanding of this regime. 
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APPENDIX 

We wish to establish that o.q* «=> 0 for n>0 is consistent with the 
Reynolds equation for q*. This result, required in arriving at Eq. (9), 
expresses the fact that the expectation of the pressure gradient may be 
written as the constant smooth value a^qg ° gi . By substitution of 
the series for q into Eq. (2), the general equation for q^^ is obtained; 

a^D^q^ o -*3a^[x"a^q^_^] -3a Vs^ “ *^n 

In terms of the Green function defined by Eqs. (10) and (11), the solution 
for is 

qi,(x) = G(x-x*)R^(x^)dx^ (A2) 


Recursive substitution of lower order q|^ on the right of Eq. (Al) making 
use of Eq. (A2) leads to a series of terms in Rp of which the general 
form is an (a - 1) - fold double integration with respect to x^,x^,..x*-, 

2 <_ a < n, of a product of (a - 1) Green functions and n factors of x* 
distriFuted between the a arguments x ,x^,...x®'. For pressure flow this 
general term may be explicitly written 



. • • |aj[x”{x®')] *'g^,|dx^dx^...dx‘'' 


times a numerical factor, depending on the set of integers p,^ which must 
a 

satisi'y 23 P ° n. For shear flow the final brace is replaced by 

(11 c: 1 

P — 1 

[x~(x*')] SpCx^ix*" )]V^/h*^| , while the numerical factor is also 
changed. 

Forming the expectation value of the integrand before performing the 
differentiation of each x factor, a valid commutation of operations, we 
obtain a CF of order n between a different sample points of the stoch- 
astic rough surface. By stationarity this CF depends on the a (a - l)/2 
argument differences. ^Carrying out next the required differentiations of the 
CF with respect to x^,x^,...x®' followed by the (a - 1) - fold 
double integration then leads to a result independent of x^, since the Green 
function derivatives also depend only on variable differences. It should be 
noted that this proof depends on the assumed rapid decay of the CF with 
distance, allowing the range of integration to be extended to infinity. 


OF PCX3R quai:ty 
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The 3. operator now remaining to the left of the expression reduces 
zero. 'Thus RJ = 0 and hence by Eq. (Al) a.q* may be chosen zero, 
choice satisfies the boundary conditions and 'establishes the desired 


it to 

This 

result. 
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